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Abstract. We recently proposed a model coupling the evolution of the 
opinions of the individual with the local network topology. The opin- 
ion dynamics is based on the Bounded Confidence model. The social 
networks is based on a group concept where each individual is totally 
connected to the members of its group and is linked to the individuals 
of the other groups with a given probability. During a time step, the 
individual has to decide between discussing with a member of its own 
network and applying the opinion dynamics, or moving groups because 
it has an opinion far from the average opinion of its own group. One of 
the main results we obtained is that the group sizes, starting from an 
homogenous situation can be strongly heterogeneous at the equilibrium 
state. This kind of heterogeneity can be identified in many real networks. 
In this paper we present the complete set of behaviours that this com- 
plex model can exhibit, at group level. In particular we will focus on the 
mechanisms that lead to the stability of the groups and the appearance 
of heterogeneity in sizes. 



1 Introduction 

Each of us knows that a system can be stable at a given scale while it is not 
at another one. That is what we are interested in and this paper is going to 
describe how population, group and individual levels can differ in stability and 
forms. But, let's come back to the beginning presenting some assumptions of the 
simple simulated model of the society we work with. 

Other's opinion is a source of cognitive inconsistency! That is what Festinger 
[B] argued adding that it is experienced as dissonance. According to him, the 
dissonance is a psychological discomfort or an aversive drive state that people are 
motivated to reduce, just as they are motivated to reduce hunger. In his balance 
theory, [TT] used a similar concept and called it imbalance. More recently, [2] 
showed that, as the dissonance and balance theories suggest, the disagreement 
from others in a group produces cognitive inconsistency and the negative states 
of dissonance or imbalance. 

The groups are a privileged place of interaction between people and the 
exchange with others can lead to dissonance. They are thus at the same time 
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the entity creating dissonance and the one reducing it. Indeed, three strategies 
can be chosen to reduce its dissonance created by the heterogeneity of the opinion 
inside its group (see |14j): changing one's own opinion to agree with others in 
the group, influencing others to change their opinions, or joining a different, 
attitudinally more congenial group. 

The two first relates to the individual interactions which are often based 
on similarity and have been extensively studied in the attraction paradigm [2] 
and other theories on interpersonal interactions as the social judgment theory 
[IB] . The third can be linked to the personal external network of the individual. 
Indeed, its external neighbours give it some information about the characteristics 
of their group and can introduce it. This stresses out the importance of the 
individual's social networks. 

The structure of social networks has been the object of study of many dif- 
ferent disciplines. In particular, a more quantitative investigation of the topo- 
logical structure of social networks has been possible with the great availability 
of data of web 2.0. The analysis performed with this tools revealed that social 
networks usually present community based structures: analyzing networks at 
different scales it is possible to identify groups of persons much more intercon- 
nected among them than with the rest of society [9 . Many different algorithm 
have been created to identify communities on large networks [S] and many mod- 
els have been proposed to explain the mechanism leading to the formation of 
such underlying structure ^5j. 

Recently many works have been done regarding evolving network topology 
and their adaptation to the social background [IQJ : as people can influence each 
other to induce a change of mind, the difference of opinion on some very impor- 
tant topics can also lead to the breaking of a social contact. In other terms, since 
people prefer to be surrounded by persons sharing similar opinion (homophily) it 
is quite likely that the change of opinions due to the opinion dynamics processes 
can lead to the change of the network structure. An interesting analysis of the 
co-evolution of opinions and networks is presented in . 

Starting from the strategies to overcome inconsistency explained above, we 
have presented in i5i a simulation model reproducing their main aspects in order 
to better understand the link between the individual choices and the organiza- 
tion of the society into groups. We will model the group concept as the one of 
community. It is based on a social network where an individual has most of its 
links to its own group and a minor part to the other groups. The interaction pro- 
cess between individuals will be modeled by one of the classical opinion dynamic 
model. 

2 The model 

The model presented in ^ couples an opinion dynamics process with an adaptive 
network structure co-evolving with the opinion. The considered network has 
a typical underlying structure based on groups. We defined group as a set of 
individuals fully interconnected among them. From the point of view of the 
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individuals, therefore, each agent has a hnk to all the other members of the 
same group. Moreover individuals can also be connected, with a fixed probability 
Pext^ to members of different communities. The resulting network therefore is 
composed by different complete graphs (the groups) connected among them 
through the links that agents establish outside the group. At the initialization 
time each agent select the membership in a group. When all the agents have 
decided their membership, the internal links are established. Secondly, each agent 
connects with probability Pext to all the agents that are not inside its own group. 
The initialization is completed randomly assigning to each agent a continuous 
opinion -d ranging in the interval [0, 1]. 

Agents choose between two possible dynamics at a time: they can commu- 
nicate with an other agent in their neighbourhood (exchanging opinion with it 
and updating their opinions with the opinion of the other agent); or they can 
reconsider their membership in the group. The choice of being part of a group, 
in fact, is usually based on homophily preferences: if the average socio-cultural 
attitudes of a group are strongly different from the one of a member, this mem- 
ber could feel a certain uneasiness and decide to change groups. In this case we 
will consider homophily preferences based on the opinion. In a simulation step 
all the agents are selected and, at their turn, they can decide with probability 
Pchange to changc their membership in a group and with probability {1-Pchange) 
to perform opinion dynamics. 

When an agent i chooses to perform opinion dynamics, it first selects a part- 
ner, j, for the contact process: the partner can be selected both among the 
members of its own group or the agents connected through external connec- 
tions. The opinion updating mechanism is the Bounded Confidence (BC) model 
introduced by Deffuant et al. in [3, and deeply analyzed in all its aspects in [1], 

The parameter e represents the tolerance of an agent and it is the funda- 
mental critical parameter of the system. The parameter /i, on the other side, 
represents the relative shift of the agents after a discussion. 

On the other side, the mechanism to change groups follows the following 
steps: the agent looks at the distance between its opinion and the average opinion 
of its own group, Oj. If \'di{t)—Oi{t)\ > e, it changes groups looking at the groups 
with which it already has a contact through an external link. It will connect to 
a chosen group with a probability: 

p ^ i-\Mt)-OAt)\ , . 

EJ^,eVi^)i^-m)-OJm ^ ^ 

where V(i) is the neighbourhood of node i. 

Once an agent change groups, its connections are completely rewired: all the 
connections he had before (inside the previous group and outside) are canceled. 
It is connected to all the agents of the new group and new external connections 
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are established with the probabihty Pext with the agents outside the new group. 
AU the relationship are symmetrical: this means that if an agent A is linked to 
an agent B, the agent B is also linked to the agent A. 

In previous papers ([4], [5]) we provided a detailed description of the opinion 
dynamics process on this kind of group-based adaptive networks, specifying the 
main differences with respect to static networks. We discovered some surprising 
behaviours compared to the known behaviour of the DefFuant bounded confi- 
dence model(BC). First of all, the coupled model exhibits a total consensus for 
an e value lower than the BC model. That is linked to the capacity of the coupled 
model to suppress the minor clusters positioned in the BC model on the extrema 
of the opinion space. In social psychology, groups are known as a source of co- 
hesion and avoidance of the isolation. Thus, that is a very interesting fact that 
the introduction of groups in the BC model suppresses the isolated individuals. 
Thus, we observed an heterogeneity of the final size distribution of groups which 
is more realistic than the more homogeneous size distribution obtained with the 
Bounded Confidence model. 

An other interesting phenomenon concernes the distribution of the group 
sizes: for certain values of £, some groups become larger as the system evolves 
while other decrease in size, sometimes until containing only one individual. 

Finally, from the point of view of the groups, the consensus remains for a 
large set of e values while, looking at the population level, there are a lot of opin- 
ion clusters. Then, each group does not simply correspond to a subpopulation 
exhibiting the same behaviour than the whole population. 

In this paper we focus on the obtained equilibrium and the form of the size 
distribution. The results presented in the following sections concern simulations 
with N — 5000 individuals and (7=500 groups. Some parameters of the model 
are fixed = 0.5, Pext = 0.001 and Pchange — 0.5) in order to study the char- 
acteristic behaviours in different ranges of the tolerance parameter e. Regarding 
e, we study the value for which the total consensus is not reached. Practically, 
it means we are interested in the various e lower than 0.267. 



3 Different tolerances, different equilibria 

Fig. [T] displays the typical evolutions of group sizes and individual opinions in 
time for one realization of various values of e. As we can observe from the right 
column, the opinion dynamics leads to the formation of stable opinion clusters 
as it happens for static networks (even if the cluster structure is not exactly the 
same as for the static case [S]). 

On the other side, looking to the left column of the plot, the situation re- 
garding groups is very different: for some values of e (0.01 and 0.06) group sizes 
continues to vary in time. Individuals, even if their opinion is stable, continue to 
jump from one group to an other and this process never reaches a steady state. 
This is not the case for e = 0.13 and e = 0.22. for which, after a transitory 
phase, groups, as individuals, become stable over time also. 
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Fig. 1. Left column: Evolution in time of group sizes. Each line represents 
one group. Right column: Evolution in time of individual opinions. Each line 
represent an individual. Each raw corresponds to a different value of e: 0.01, 
0.06,0.13,0.22 
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We can argue therefore the presence of different kinds of equihbrium, for 
some vakies of the tolerance parameter e, we find stable opinions and unstable 
groups, above a threshold value, also the groups stabilize. To study this transition 
to static equilibrium we consider the relative standard deviation of the size in 
the T=100 last time steps at the population level opinion equilibrium, for each 
of the groups: 

{Si)t 

where the symbol (...)t indicates the average on time. 

This indicator, averaged for all the groups and on 100 independent realiza- 
tions is displayed, with its relative error, in Fig. [5] As we can observe from this 
figure, the transition to group stability can be set at £ ~ 0.1. 




E 



Fig. 2. Relative standard deviation on time, averaged on all the groups and on 
100 realizations 

This result is coherent with what we observed in [5] where we studied the 
percentage of individuals susceptible to change groups as a function of e. For 
e < 0.1 it is impossible to obtain the situation one opinion cluster in one group; 
therefore the dynamics of groups cannot converge to a steady state. This situa- 
tion is described in Fig. [31 First of all, we can notice that the number of opinion 
clusters inside a group is always lower than the number of clusters at population 
level. In particular, for e ~ 0.1, we observe that the average number of opinion 
clusters at population level is five, but each group contains only one cluster. 
For £ < 0.1, therefore, the system is unstable: when there are several opinion 
clusters by group, some individuals are necessary far from the average opinion 
of the group since the distance between two opinion clusters is equal to at least 

£. 
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Fig. 3. Number of opinion clusters by group and number of opinion cluster at the 
population level for various values of e. The result is averaged on 30 realizations. 
The dotted line represents the minimum and the maximum on the realizations. 



Since in the following we are going to explore some statistical properties of the 
size distributions, it is important to analyze in details the non-equilibrium phase. 
As we observed before, for e < 0.1 the groups are not stable. Notwithstanding 
this dynamical instability, if we observe the aggregate properties of the system, 
namely the size distribution at different times of the evolution, we can argue 
that the statistical properties of the model remain constant. In the Fig. H] we 
analyzed the case e — 0.01 and we displayed the group size distribution at 
different moments of the evolution. The figure shows that the size distributions 
maintain over time the same shape. This confirms that the size distribution 
remains a good indicator of the heterogeneity of the group size in the non- 
equilibrium phase. 



t=500 



t=1000 



t=1500 




Fig. 4. Group size distribution for e = 0.01 at different moment (T=500, 
1000,1500) of the evolution. The result is averaged on 100 realizations. 
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4 The various levels of heterogeneity 

From the Fig. [T] we can also extract interesting considerations regarding the 
distribution of group sizes. Since we use a random procedure to fill the group 
structures, at the initial state, all the groups are distributed with a small devia- 
tion {as = 3.15) around the average value: (S) = N/G = 10. The average size of 
the larger group, at the beginning is Smax = 20.8 ± 1.62. Of course the average 
size of the groups will remain the same during all the evolution of the system 
since the population N and the number of groups G are constant. By the way, 
looking at the Fig. [T] and taking into account the considerations about Fig. HI 
we observe that for e = 0.01, some bigger groups {S ^ 100) already appear. At 
e = 0.06 a macroscopic group containing more than 20% is present, while all 
the other groups seem to be much smaller. The situation is again different at 
£ = 0.1 where many intermediate sizes seem to be allowed. Finally, for e ~ 0.22 
the situation seems to be similar to the initialization moment. A statistical de- 
scription of these different situations is displayed in Fig. [S] where the cumulative 
size distributions for the four cases analyzed in Fig. [T] are reported. 




Fig. 5. Cumulative size distribution at the end of the evolution for e = 
0.01,0.06,0.13,0.22. The dotted lines represent the size distribution at the ini- 
tialization. The result is averaged on 100 realizations. 

Fig. [S] provides a statistical description of the assumption we did by the 
observation of the typical run in Fig. [T] for e — 0.01 group sizes present a 
distribution that decrease fast with the size and that can be fitted with an 
Exponential law. For e = 0.06 we observe a bi-modal distribution, where the 
large tail corresponds to the few macroscopic groups. For e = 0.13 group sizes 



How opinion dynamics generates group hierarchies 



9 



present a sort of hierarchy structure where ah the sizes seem to be possible. 
The data, in this case, can be fitted with a power law distribution, denoting the 
absence of a typical scale in the system. Finally the size distribution for e ~ 0.22 
is not strongly dissimilar to the initial one. 

A more complete set of cumulative degree distributions for different values 
of £ is given in figure |6l 




10 100 JOOO iOOOOJO JOO 1000 10000 iO 100 1000 10000 10 100 1000 10000 10 100 1000 lOOOO 



Fig. 6. Cumulative size distribution at the end of the evolution for all the ana- 
lyzed values of e. The result is averaged on 100 realizations. 



The variability of the sizes, at the final state is described from the second 
moment of the size distribution {S'^)i,reai- We can provide a characterization 
of the heterogeneity level of the group structures introducing the parameter: 
K — (5'^)/((5))^: if K 1 the variance is dominated by the second moment 
and therefore the sizes are strongly heterogeneous, while for k ^ 1 we have an 
homogeneous distribution. The Fig. [7] displayed this parameter as function of e. 
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Fig. 7. Heterogeneity between group sizes as a function of e. The result is aver- 
aged on 100 realizations. 

As we can see, the parameter n quickly increases for e > 0.02, reaches a pick 
for e ~ 0.07 and rapidly decreases. Connecting this information with the analysis 
of the size distributions (Fig. ^ we can argue that: for e < 0.02 the system 
presents an exponential size-distribution and, therefore, it is characterized by 
a low value of heterogeneity. The range 0.02 < £ < 0.1 corresponds to the bi- 
modal distribution, with the two characteristics scales (small groups and few 
macroscopic groups). This is the phase of higher heterogeneity. The range 0.1 < 
£ < 0.14, instead, is characterized by a scale free group size distribution. Also this 
phase presents a high level of heterogeneity, even if lower than in the previous 
case. Finally the distribution returns gradually to the initial one (passing through 
situations very similar to the ones obtained for e < 0.02). For £ > 0.2 the groups 
result to be homogeneous in size among them. 

5 An heuristic explanation of stability and heterogeneity 

The mechanism that generates heterogeneity is a hidden preferential attachment 
that estabilishes the rapid growth of some particular groups. The model in fact 
implies that an individual, when he chooses a group, considers in his selection 
only the groups on which he can retrieve information, namely the groups in which 
he has a friend (a link). Since the choice of the external links does not depend on 
the membership (an agent has the same probability pext to select all the agents 
that are non member of its own group), the probability that an external link 
ends in a big group is higher. The fact that more links bring to the bigger groups, 
makes these more eligible as possible destinations when agents decide to change 
groups. The creation of heterogeneity is therefore an effect of this preferential 
attachment. The different shapes of the size distribution (exponential, bi-modal, 
power law) are a consequence of the fact that for different values of e the rich 
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get richer mechanism happens in different ways. At the beginning, due to the 
random initiahzation, the average opinion of each group is 0/ 0.5. Therefore, 
on average, ah the individuals with an opinion z9 situated in the range 0.5 — e < 
1? < 0.5 + e are stable: since their opinion is very similar to the average opinion of 
the group, they are satisfied and they won't change groups. The group dynamics 
therefore starts from the agents with opinions further than e from the center 
of the opinion interval. When e — the number of agent that are initially 
susceptible to move is higher. 

Another important point is that the average initial number of individuals 
in a group is: ni — N/G. The average distance between two individuals at the 
initial point therefore is: d = l/{ni — 1), since the opinion space measures 1. 
When d > e opinion dynamics, on average, does not happen. With our choice of 
parameters it is always the case for e < 0.1. 

Two mechanisms linked to the both dynamics are responsible for the hetero- 
geneity emergence. At a first stage, the initial small unbalance due to stochas- 
ticity of the average opinion of groups, favors a polarization of the groups. The 
Fig [5] shows how it occurs due to the dynamics of group change. It represents 3 
groups with 6 agents (grey squares) in each. Agents with an opinion far from the 
average (the dark vertical line) start to move to groups with a closer O/. These 
moves implies a polarisation of the average opinion of the group which, at its 
turn, leads in each group to a desertification of one opinion extreme in favour of 
the other. 
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Fig. 8. Group dynamics polarisation effect. Evolution of a population of 18 in- 
dividuals in 3 groups. The time goes from the left to the right starting from the 
top. Time to 3. 



As people condense on one extreme, opinion dynamics happens reducing 
the standard deviation of the group opinion. The Fig. [S] shows how two close 
(regarding e) individuals infiuence each other leading to more cohesiveness of 
the group members because the minimum opinion of the group become closer to 
the maximum one. This standard deviation decreasing phenomenon gives to the 
new individual arriving in the group more chance to be satisfied and to remain 
in the group. 
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Fig. 9. Opinion dynamics decreasing standard deviation effect. Evolution of a 
population of 18 individuals in 3 groups. The time goes from the left to the right. 
Time 4 and 5 



One can see that at a given time, the number of individuals in each group can 
differ a lot. That is sufficient for the preferential attachment described previously 
to allow a larger group to become larger and larger. Depending on the number of 
opinion clusters, the more or less remaining central people are going to attract 
the more extreme opinion clusters and, if they are enough, make the group 
having an average opinion close to the center. That is the most frequent case. 
The exception is when e leads to three clusters at the population level. 

This phenomenon leading to the emergence of heterogeneity does not occur 
completely for e — because all the agents (also the central ones) move since 
the beginning, forbidding the polarization of group opinions. For higher values 
of e the polarization starts to occur, but the probability of finding a group 
able to gather the central opinions ( i.e. a group containing a higher number of 
central individual stable since the beginning) is lower, since the stability interval 
is very small. Just one or few groups have this characteristic and, therefore, the 
capacity of starting the rich get richer mechanism. This leads to the bi-modal 
distribution. 

Finally, as e increases, the number of groups suitable for the preferential 
attachment also increases, allowing the simultaneous growth of more than one 
macroscopic group. This is the origin of the power law distribution. 

Stability is connected to the relative speed of the opinion dynamics process 
and of the group dynamics. As e — > opinion dynamics became slower, never 
allowing the formation of one cluster in one group. Opinion dynamics, in fact, for 
low values of e can happen only at population level on large time scale, without 
influencing the group formation. As we observed before, for higher values of e 
opinion dynamics start to be possible also at group level. Since the number of 
connections inside groups is higher, this mechanism accelerates opinion dynamics 
process. The combination of opinion dynamics and group dynamics bring in 
such cases to specification of opinions inside groups, one cluster in one group. 
In intermediate cases, like for e = 0.06, many small groups composed of two 
symmetrical side clusters remains. This situation allows a continuous exchange 
between the central cluster and the groups in such situation. 
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6 Conclusions 

In this paper we analyzed the properties of the groups obtained by a combined 
dynamics of opinions and membership choices. In particular we analyzed two 
aspects: the different kind of equilibrium and the size distributions. 

For what concerns the equilibrium, two main situations can be identified: for 
e < 0.1, even if all the opinion clusters are formed, the dynamics of groups never 
reaches a stable state. In this case, however, even if the individuals continue to 
jump from one group to the other, at the aggregate level the degree distribution 
results stable. For e > 0.1, instead, the system reaches an equilibrium where 
neither the opinions nor the groups change over time. 

Regarding the group sizes we identified four characteristic regions: low het- 
erogeneity phase (e < 0.02); few groups dominance (0.02 < e < 0.1) where 
macroscopic structures are formed; group hierarchy (0.1 < e < 0.2) and group 
homogeneity (e > 0.2). 

The phenomenological analysis we provided, by the way, highlights some 
important points: first of all, the dynamics strongly depends on the number of 
initial groups and, in particular from the initial number of individuals in the 
groups. This point is a subtle one and a more analytical study would be needed, 
through a mean field analysis. 

Moreover, we tried to eliminate the preferential attachment mechanism, giv- 
ing to the agents the possibility to do a rational choice between all the groups 
independently from the network structure. With this option, we didn't observe 
any heterogeneity for £ < 0.1 but some heterogeneity remains in the other case. 
Also in this case, some simple mathematical models, would be needed to provide 
a better quantitative explanation. 

Finally, what we obtained, except for very low values of e, is that the het- 
erogeneity level of sizes decreases when e increases. The significance of such a 
result is not easy to catch. Indeed one can notice that the e used for group as- 
sessment and the one used for individual assessment are quite different in their 
meaning. They could be considered as two parameters in another model. On 
the one hand, smaller is e, larger is the need for cohesion of the individual to 
recognize a group as its own group. On the other hand, larger is e, more tolerant 
is the individual regarding the opinion difference to others. Concretely, it means 
that an individual more defines itself by its own opinion when e is low while it 
more defines itself by its group membership when e is large. One can see a kind 
of individualistic society versus a collective one. That is very counterintuitive is 
that the individualistic society leads to the emergence of large majority groups 
while the collective society leads to a much more homogeneous size of groups! 
A deeper investigation in the social psychological litterature would be necessary 
to assess about the realistic feature of this result. 
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